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Abstract: We compute the one-loop effective potential for noncommutative U(l) gauge 
fields on S^xS^. We show the existence of a novel phase transition in the model from the 
4— dimensional space x S| to a matrix phase where the spheres collapse under the effect 
of quantum fluctuations. It is also shown that the transition to the matrix phase occurs at 
infinite value of the gauge coupling constant when the mass of the two normal components 
of the gauge field on S? x S? is sent to infinity. 
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1. Introduction 

Fuzzy approximations of spacetime (like lattice regularizations) are designed for the 
study of gauge theories in the nonperturbative regime using Monte-Carlo simulations. They 
consist in replacing continuous manifolds by matrix algebras. The resulting field theory 
will thus only have a finite number of degrees of freedom and hence it is regularized. The 
claim is that this method has the advantage -in contrast with lattice- of preserving all 
continuous symmetries of the original action at the classical level || ||. 

Field theory on the fuzzy sphere is the most studied example in the literature. In 
perturbation theory it is shown that scalar field theories on suffer from the UV-IR mixing 
problem || . Moreover it is shown that there exists new nonperturbative phenomena which 
are missing in the commutative theory. For example a novel phase has been discovered 
in scalar field theories on (the so-called non-uniform phase or matrix phase) which 
has no commutative analogue Ml. This new phase was also observed in three dimensions 



|1C]. Generalization to 4— dimensional fuzzy spaces and their scalar field theories were 



undertaken in [11]. 

The quantum properties of the gauge field on the fuzzy sphere have been studied in 
0, |, |. In Q the effective action was computed to one loop for U(l) gauge fields. It was 
shown that the model contains a gauge invariant UV-IR mixing in the limit L — > 00, i.e 
the effective action does not go over to the commutative action in the continuum limit. 
Furthermore a first order phase transition was observed at one-loop from the fuzzy sphere 
phase to a matrix phase where the sphere collapses. This transition was previously detected 
in Monte Carlo simulation of the model reported in ||. In some sense the one-loop result 
for the U(l) model is exact. 
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It was also shown in []J that the UV-IR mixing and the matrix phase both disappear 
in the limit where we send the mass of the normal scalar component of the gauge field 
on S| to infinity. This means in particular that the nonperturbative S^-to-matrix phase 
transition is a reflection of the UV-IR mixing seen in perturbation theory and that this 
latter finds its origin in the coupling of the normal scalar field to the two dimensional gauge 
field. The differential calculus on the fuzzy sphere is intrinsically 3— dimensional and as a 
consequence there is no a gauge-covariant splitting of the 3— dimensional fuzzy gauge field 
into its normal and tangent components on S^; hence the action will necessarily involve 
the interaction of the two fields. This result (among many others) was confirmed recently 
in our Monte Carlo simulation of the model where we have also found a novel third-order 
one-plaquette-like phase transition which the model undergoes and which we can also trace 
to the coupling of the normal scalar field. The full phase diagram of the model will be 
reported elsewhere 

The main goal of this article is to study the phase structure of £7(1) gauge theories on 
fuzzy S^xS^. The advantage of considering S? x S? is that one can use all the machinery 
of the well known SU (2) polarization tensors. Other studies of noncommutative gauge 
theories on 4— dimensional fuzzy spaces have already appeared [12|. 

This article is organized as follows. In section 2 we give a brief description of the 
geometry of fuzzy S^xS^. In section 3 we introduce fuzzy gauge fields and we write down 
the action we will study in this article. In section 4 we compute the effective potential. In 
section 5 we show the existence of a first order x S^-to-matrix phase transition in exact 
analogy with the two-dimensional case and we derive the critical line. Section 6 contains 
the conclusion. 



2. Fuzzy S| x S| 

Fuzzy x S| is the simpliest 4 dimensional fuzzy space. It is a finite dimensional 
matrix approximation of the cartesian product of two continuous spheres. This fuzzy space 
is defined by a sequence of Connes triples |j| 

S 2 L xS 2 L = {Mat {L+1) 2,H L ,A L }. (2.1) 

Mat^ L+1 ^2 is the matrix algebra of dimension (L+ l) 2 and is a suitable Laplacian 
acting on matrices which encodes the geometry of the space. It is defined by 

A L = [L AB ,[L AB ,-]]=C\ B (2.2) 

where Lab, with A, B = 1, 4, are the generators of the irreducible representation (-|, |f) of 
SO (4). The generators Lab (with Lab = —Lb a) satisfy the commutation relations 

[L ab ,Lcd] = JabcdefLef 

= $bc'Lad — $bdLac + &adLbc — Sac'Lbd- (2.3) 

Hl in (|2.1| ) is the Hilbert space (with inner product < M,N >= jp-^x Tr (M^N)) which 
is associated with the irreducible representation (-§,-§) of SO (4). 
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Since 50(4) = [5*7(2) x 5*7(2)] /Z 2 we can introduce 5*7(2) (mutually commuting) 
generators L)p and by -2Ly = \e a bc^bc + L ai and -2L^ = \e ahc L hc - L a4 with 
a = 1, 2, 3 and e a b c is the three dimensional Levi-Civita tensor. Then it can be easily shown 
that the two 50(4) quadratic Casimir can be rewritten in the form (where cabcd is the 
four dimensional Levi-Civita tensor) 



T 2 

L AB 

cabcdLabLcd 



4[(4 1 )) 2 + (4 2 )) 2 ]=2L(L + 2)^8c 2 

W$>? - (4 2) ) 2 ]=o. 



(2.4) 



Similarly the Laplacian C\b rea ds in terms of the three dimensional indices as follows 



r 2 

L AB 



Hi) 



+ 



(2.5) 



where Ch = [ii 1 ', • ] and = 
by the coordinate operators 



L 1 ^',- ]• For S£ x S£ the algebra Mat^ L+1 ^2 is generated 



which satisfy 



R 



r(l) 



R 



r(2) 



(2.6) 



E 

a=l 



-(0 



R?l, 



,0) ~0') 



-Oij£abc% c i 



1,2. 



(2.7) 



In the limit L — > oo keeping i?i and i?2 fixed we recover the commutative algebra of 
functions on S 2 x S 2 . If we also choose to scale the radii R\ and R 2 such as for example 
0\ = R\/L\ and # 2 = R\jL 2 are kept fixed we obtain the non-commutative Moyal-Weyl 
space M.q 
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The algebra of matrices Mat( i+1 p can be decomposed under the action of the two 
SU(2) of 50(4) as Mat^+i) ® Matf L+ iy As a consequence a general function on S 2 x S| 
can be expanded in terms of polarization tensors |l4| as follows 



^kxmxk^m-iXkx 

fci=0mi=— fei ^2=0 m,2=— A;2 



E 



mi 



3. Fuzzy gauge fields 

U(n) gauge field on S|xS 2 can be associated with a set of six hermitian matrices 
Dab £ ^ a t n {L+i) 2 (Dab = —Dba) which transform homogeneously under the action of 
the group, i.e 

Dab - UD AB U-\ U e U {n(L + l) 2 ) . (3.1) 

In this paper we will be mainly interested in *7(1) theory on 5£ x 5 2 . The action is given 
by (with Ttl = 1 . 2 Tr, g is the gauge coupling constant and m is the mass of the normal 

(L+l) 

components of the gauge field ) 
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S - ^ -7^1 [Alii- A'/)]" + - 1 ALU' DLi--'I'r L \D. \1J-D ( ■!>][) h:i .- 



2 2 

+ ^T^ Tr dD\ B ~ L\ B ? + ™ 2 Tr L {e ABC DD AB D CD f. (3.2) 



The equations of motion are given by 

i[D C D,F AB ,cD} + ^{D AB , $i + $ 2 } + ^{cabcdDcd, $1 - $2} = 0. (3.3) 

VC2 v/C2 



As we will see shortly F ABj qd = i [Dab, F>cd] + IabcdefDef can be interpreted as the 
curvature of the gauge field on fuzzy S| x whereas $1 and $2 (defined by D AB — L AB = 
8^/02(^1 + $2) and € AB cdD ab Dcd = 16y / C2(^ ) i — $2)) can be interpreted as the normal 
components of the gauge field on S| x S^. 

The most obvious non-trivial solution of the equations of motion (^^) must satisfy 
F AB ,cd = 0, D AB = L 2 AB and € ABC dD ab D C d = (or equivalently F AB = 0, = 0). 
This solution is clearly given by the generators L AB of the irreducible representation (■§,■§) 
of SO (4), viz 

D AB = L AB . (3.4) 

As it turns out this is also the absolute minimum of the model. By expanding D AB around 
this vacuum as D AB = L AB + A AB and substituting back into the action ( [3 .2D we obtain a 
U(l) gauge field A AB on x S| with the correct transformation law under the action of the 
group, namely A AB — >U A AB U~ l + U C AB U~ l . The matrices D AB are thus the covariant 
derivatives on S^xS|. The curvature F A b,cd i n terms of A AB takes the usual form 
F ab ,cd = iC-abAcd ~ iC-CDA AB + f ABC defA E f + i[A AB , A C d\- The normal scalar fields 
in terms of A AB are on the other hand given by 8-^/02 (3>i+<I>2) = L AB A AB + A AB L AB + A 2 AB 
and 16 % /ci($i - $2) = (-ABCd{LabA C d + A ab L C d + A ab A C d)- 

We can verify this conclusion explicitly by introducing the matrices D a ^ = L a + A a 
and D { a ] = L { a ] + A a 2) defined by 



D^= - - 

a 2 



1 

-j£abcF>bc + F>oA 



a - 2 



1 

-jtabcF>bc - D a 4 



(3.5) 

Clearly Djp O^a ) and .D^ (-A^) are the components of D AB (A AB ) on the two spheres 
respectively. The curvature becomes F^ = iC$A^ -iC^A$ + (%e a6c A£* ) + i[A^, A^} 
whereas the normal scalar fields become 2 A /c 2 "<I>j = (D^) 2 -c 2 = L$A$+A$L%' + (A%') 2 . 
In terms of this three dimensional notation the action (|3.2j) reads 



S = S^+S^ + ^- 2 Tr L {F^) 2 . (3.6) 

and are the actions for the U(l) gauge fields A a 1 ^ and A^ on a single fuzzy sphere 
S^. They are given by 



^ir ) 4' , -^[4• | ■4' , ]4 i, 



2m 



2 



+ — Tr^<I> 



2 



(3.7) 
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It is immediately clear that in the continuum limit L — >oo the action ( |3,6| ) describes 
the interaction of a genuine 4— d gauge field with the normal scalar fields $j = ni A$ 

(i) 

where n a is the unit normal vector to the i-th sphere. The parameter m is precisely the 
mass of these scalar fields. Let us also remark that in this limit the 3— dimensional fields 
Af decompose as A® = (A^) T + r$$i where (A$) T are the tangent 2— dimensional 
gauge fields. Since the differential calculus on x S| is intrinsically 6— dimensional we 
can not decompose the fuzzy gauge field in a similar (gauge-covariant) fashion and as a 
consequence we can not write an action on the fuzzy x S| which will only involve the 
desired 4— dimensional gauge field. 



4. Quantum effective potential 

The partition function of the theory depends on 3 parameters, the Yang-Mills coupling 
constant g, the mass m of the normal scalar fields, and the size L of the matrices, viz 

r 4 

Z L [J,g,m]= H [dX AB ]e~ s ^ +Tr ^ x ^. (4.1) 
■> A<B=1 

In the background field method the field is decomposed as Xab = Dab + Qab where Dab 
is the background we are interested in studying and Qab stands for the fluctuation field. 
We add the usual gauge fixing and Faddeev-Popov terms given by 



S g . f + S gh = -J-^Tr L [DAB 'Q AB]2 + ^TtlcIDab, [D AB ,b]]. (4.2) 



Performing the Gaussian path integral we obtain the one-loop effective action 

T [D AB ] = S [Dab] + ^Tr e TRlogQ A BCD - TR\ogV 2 AB . (4.3) 



Qabcd is defined by 



L ->2 _ x (-, _ lN \ -n , -n o*t,„„„ _l 4m nC 1 ) 



^ABCD = -^EF^ABflD ~ \1 ~ | J T^AB^CD ~ 2iTaBCD + ^2 ABO D ' ( 4 - 4 ) 

where 5ab,cd = ^ac^bd - SadSbc, and 

^abcd = (D EF — L ef )5ab,cd + -^{^efghDefDgh)^abcd 

-VabVcd - VabVcd + ADabDcd + ±D ab D C d. (4.5) 

The notation T>ab and Tabcd means that the covariant derivative Dab and the curvature 
Fabcd act by commutators, i.e V A b(M) = [D A b, M], T A bcd{M) = [F ABC d , M] where M 
is an element of Mat^ L+ ^2. Wwe have also introduced the notation Dab = \^abcdDqd- 
TR is the trace over the 4 indices corresponding to the left and right actions of operators 
on matrices. Tr§ is the trace associated with the action of SU(2)xSU(2). 

The main goal of this article is to check the stability of the solution ( |3.4| ), in other 
words to check whether or not the fuzzy space x is stable under quantum fluctuations. 
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Towards this end it is sufficient to consider only the background field Dab = ^Lab where 
the order parameter <p plays the role of the radius of the two spheres of S^xS^. Therefore 
the computation of the effective action reduces to the computation of the effective potential 
V e s((f)) = T[cPLab}- The classical potential is given by 



(4.6) 



V^Lab] = - \? + \m 2 (<? - I) 2 ) 

The effective potential (in the gauge £ = 1) is given by 

V eB = V + ^Tr 6 TRlog<p 2 - TRlogtf + ^Tr$TR log &abcd 

= V + 4(L + l) 4 log</> + ^Tr 6 TRlo g n A BCD. (4.7) 
We are only interested in the (^—dependence of the operator Q which is defined by 

&ABCD = -^C 2 ef 8ab,CD + — — ^ f ABCDEF&EF + jp. — ^ABCD' ( 4 - 8 ) 

where 

^■ABCD = f 1 ~~ L 2 ef 5aB,CD — C-AB&CD ~ ^AB^CD + ^AB^-CD + ^AB^CD- 

(4.9) 

We will need to use the following identities 

XabYab = 4(x( 1 )y« + x^yP) , 

fABCDEFTr [XabYcdZef] = We abc Tr [x^Y b {l) + X^Y b {2) ' Z&\ . (4.10) 

The matrices Xa (la ) are related to the matrices Xab (Yab) by equations of the form 
Q3.5| ). Using these identities we can express the last term in ( f4.7D in the following way 

^Tr 6 TRlo g n ABCD = J dXABe- 7 ^^ 1300 * 00 

dX^ e -*TrxPn ab xV 

= Tr 3 TR log n ab . (4.11) 

The contributions coming from the two spheres are equal and hence the factor of 1 (instead 
of I) in front of the last logarithm. Tr% is the trace associated with the action of SU(2) 
on the two dimensional sphere. The Laplacian £l ab is defined by 

n ab = 2C 2 AB 5 ab + 16 (l - ij ie abc 6p + 8m 2 n^ b \ 
fi2=4Pi 1) -l4 1) 4 1, +2(l-^)^. (4-12) 
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pfv* is the normal projector on the fuzzy sphere defined by P 



(i) 



ab 



x^x^ where are 



the coordinate operators defined in ( |2.6| ) with R\ = R2 = 1. The presence of this projector 
means in particular that we can not diagonalize in the polarization tensors basis. However, 
in order to have an idea of the phase structure of the model, we can expand around m = 0. 
This approximation was more than sufficient in the two-dimensional case as discussed in 
great detail in fij. Therefore it is convenient to separate the logarithm term as 



logO Q 



logO^ + log 



1 + 8m 2 



1 



o 



(4.13) 



O^? is clearly equal to £l ab when m? = 0. This operator can be trivially diagonalized 
in the vector polarization tensors basis (Y^ 1 1 ) a on the first sphere tensor product the 
scalar polarization tensors basis Yi 2m2 on the second sphere. Indeed by introducing the 
total angular momentum on the two-dimensional sphere = Ca 
the generators of SU{2) in the spin 1 irreducible representation we can rewrite O ab 
following form 



^ where are 



>(0) 



in the 



-o (0) 



(£P) 2 5 ab + (£?>f5 ab -(l 



(2)\2. 



[{jP) 2 ab -{C^Y5 ab -25 ah ] 



Hence it is convenient to use the following expansion for the matrices in ( 4. lip 

*P= E Y,€T(nr) a ®^. (4.i5) 

Thus 



(4.14) 



Tr 3 TRlog 0$ = Y, (2ji + 1) (^2 + 1) log 



1 



+ 1 ) + . 



2 1^2 



+ 1 



(4.16) 



In the limit L — >oo it is easily verifiable (for example by making an expansion in 1 — ^) 
that this term is subleading compared to L 4 . The second contribution in the limit m — >0 
is given by 

1 \ ~~ o_ ( 1 



Tr 3 TRlog ( l + 8m* 
8m 2 



C2 



O(o) 
Tr 3 T.R 



(1) 



1 



O(0) 



32m z Tr 3 TR 



dp6 b l) + 16m^ I ! 



O(0) 



(l)J 1 ) 



1 



Tr 3 TR 



1 



(4.17) 



In the large L limit it is possible to show (see the appendix) that all terms in ( (4.17|) are 
subleading compared to the L 4 behaviour seen in the second term in ( |4.7|) and hence the 
full one-loop quantum contribution to the effective potential is given by the logarithmic 



potential in (4.7). Thus as long as we are in the region of the phase space near m ~ the 
effective potential behaves in the large L limit as follows 



1 



VeS _ 

4L 4 Ag 2 L 2 \A 



1 JA - V + ± m 2 ,.2 



1 

—'l 

4 



1) 2 )+Iog0. (4.18) 
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This result is to be compared with the quantum effective potential for U(l) gauge fields 
on a single fuzzy sphere S 2 computed in GJ which is given explicitly by 



VeS J_ 

L 2 2g 2 



1 



+ i m 2 (A? 



1) 



+ log< 



(4.19) 



V eS (<P) = 2c 2 A^ 2 a 4 



+ 4c2 log (j) + subleading terms. 



(4.20) 



5. The S|x S|-to-matrix phase transition 



The second term in the potential ( 4.18 ) is not convex. This implies that there is a 
competition between the classical potential and the logarithmic term which depends on the 



values of m and g. The equation of motion ^j^- = will admit in general two real solutions 
where the one with the least energy can be identified with the fuzzy S^xS 2 solution (|3.4|), 
This equation of motion reads 



(1 + m 2 



m 2 4) 2 + ig z L 



2 T 2 



0. 



The quantum solution is found to be very close to 1, viz 

4g 2 L 2 



1 



1 + 2m 2 



+ 0((g 2 L 



2 T 2\2\ 



(5.1) 



(5.2) 



However this is only true up to an upper value of the gauge coupling constant g (for every 
fixed value of m) beyond which the equation of motion ceases to have any real solutions. 
At this value the fuzzy S 2 xS 2 collapses under the effect of quantum fluctuations and we 
cross to a pure matrix phase. In other words we can not define a gauge theory everywhere 
in the phase space. As we will see below when the mass m is sent to infinity it is more 
difficult to reach the matrix phase and hence the presence of the mass makes the fuzzy 
Sf x S| solution ( fOj) more stable. 

The critical value can be computed by requiring that both the first and the second 
derivatives of the potential V e R with respect to cj> vanish. In other words, for every fixed 
value of m the critical point is defined at the point (g*,m) of the phase space where we go 
from a bounded potential to an unbounded potential. Solving for the critical value we get 
the results 



8(1 + m 2 



l + \ 1 + 



32m 2 (l + m 2 



and 



29tL 2 = - 1 -(l + m 2 )<Pi + 1 -<Pl + ^<Pi 
In the particular case of m 2 = the critical value is 



2 T 2 



9tL 



1 / 3 



2 V8 



(5.3) 



(5.4) 



(5.5) 
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Extrapolating to large values of the mass (m — >oo) we obtain the scaling behaviour 



a 2 T 2 



m 



16 



(5.6) 



In figure 1 we plot the phase diagram defined by this equation 1 . As we increase the value 
of the coupling constant g (for a fixed value of m 2 ) there exists a critical point where the 



fuzzy x solution becomes unstable and thus the minimum (3.4) disappears. Similarly 
as the value of the mass squared m 2 increases (for a fixed value of the coupling constant 
g) there is a critical point m 2 where x collapses. Clearly the value of m 2 is found 
by inverting equation (|5.6|) , viz 



l6g 2 L 2 + 1-V2. 



(5.7) 



Finally we remark that as the value of m 2 increases it is more difficult to reach the transition 
point, in fact when m 2 — >oo the critical value g 2 approaches infinity. 



The phase diagram 







equation (5.6) + 




matrix phase 








fuzzy sphere phase 









0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.1 



Figure 1: The S^xS^-to-matrix critical line. 



6. Conclusion 

We have described the qualitative behaviour of a first order phase transition which 



occurs in a U(l) gauge theory on x S 2 L . Using the one-loop effective potential ( 4.1 g ) 
of this theory we found that there exists values of the gauge coupling constant g and the 
mass m for which the fuzzy x S| solution (|3.4| ) is not stable. Thus for these values a 
U(l) gauge theory on x is not well defined. This means in particular that the model 



1 Notice that if we allow m 2 to take negative values, the gauge coupling constant g 2 will be a more 
complicated function of m 2 . However we are only interested in positive values of m 2 for which the behaviour 



(5.3) and (5.4) 



of <?» as a function of m is the straight line (5.6) which can be deduced from the large m behaviour of 
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( |3.2| ) can be used to approximate U(l) gauge field theories on S 2 x S 2 only deep inside the 



fuzzy sphere phase. However it is obvious from the critical line (5.6) that when the mass 



m of the two normal scalar fields on S 2 ^ x S 2 goes to infinity it is more difficult to reach 
the transition line. Therefore we can say that our main goal of defining a nonperturbative 
regularization of a U(l) gauge theory on S 2 x S 2 is achieved. Generalization to U(n) with 
and without fermions should be straightforward as long as we are only interested in the 
effective potential. 

Acknowledgements The authors P. Castro- Villarreal and R. Delgadillo-Blando would 
like to thank Denjoe O'Connor for his supervision through the course of this study. B. Ydri 
would like to thank Denjoe O'Connor for his extensive discussions and critical comments 
while this research was in progress. The work of P. C. V. was supported by CONACYT 
Mexico Grant No.43891-Y. The work of R. D. B. is supported by CONACYT Mexico. 

A. Evaluation of 1%, I 2 and J 3 . 



In this appendix we show that the 3 terms in ( 4.17| ) are subleading compared to L 4 . 
Let us define 

h = ±Tr 3 T R (J-) , h = ±Tr 3 T R (J-) 



We evaluate these traces by using the base of polarization tensors ( Y^ 1 1 J ®i^ 2m2 . Usin 

)ab 



the identity C a Cb = C 2 5 a i, — (8 ■ C) a b — (9 ■ C) 2 ab the eigenvalues of the operator C a Cf, are 



given by 



m xh = \ tii tii + 1) - h (h + !)) 2 " \ Ui til + + k (h + 1)) , (A.2) 

5(0) 



whereas the eigenvalues of Q ab given by ( 4.14| ) are 



4T = 8 (*i (4 + 1) + ^2 (4 + 1)) + 8^ til til (h + 1) - 2) • (A.3) 

Using these facts the two quantities I\ and I 2 can be shown to be given by 



1 



In order to evaluate I3 we notice the fact that is proportional to ^Y 00 ^ thus by using 
the algebra of vectorial polarization tensors we get the identity 

Tr { (Y\ M ■ Y°°) (Y?° • Y+ } = (L + l)(2i + 1) j \ | { \' . (A.5) 
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The final result for ^3 is 



2 



2U + 1 



X 



(L + 1) (2*i + 1) 













h 


j] 






L 




1'] 


I* 


2 







(A.6) 



The large L behaviour of I\ , I2 and I3 can be studied with the help of the different identities 
of [14]. The first sum in ( |A.4| ) diverges at most as I? in the continuum L — > 00 limit and 



hence I± converges to zero as 1/L 2 . On the other hand the sum in I2 behaves at most as 
L 4 thus the whole expression goes to zero as 1/L 2 . For J3 we can check that the sum goes 
as L, i.e I3 appraoches as 1/L 3 . 
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